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Abstract. This paper investigates tlie existence conditions of cusp points in the design parameter 
space of the RPR-2PRR parallel manipulators. Cusp points make possible non-singular assembly- 
mode changing motion, which can possibly increase the size of the aspect, i.e. the maximum sin- 
gularity free workspace. The method used is based on the notion of discriminant varieties and 
Cylindrical Algebraic Decomposition, and resorts to Grobner bases for the solutions of systems of 
equations. 
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1 Introduction 

It is well known that the workspace of a parallel manipulator is divided into 
singularity-free connected regions (2). These regions are separated by the so-called 
parallel singular configurations, where the manipulator loses its stiffness and gets 
out of control. The so-called cuspidal manipulators have the ability to change their 
assembly-mode without running into a singularity, which thus may increase the size 
of the singularity-free regions 19] |3- The word "cuspidal" stems from the notion 
of cuspidal configuration, defined as one configuration where three direct kinematic 
solutions coalesce. A cuspidal configuration in the manipulator joint space allows 
non-singular assembly-mode changing motions. Thus, determining cuspidal config- 
urations is an important issue that has attracted the attention of several researchers 
m |6l [T3] [T] . In particular, llT3l (resp. HI) has analyzed the cuspidal configurations 
of planar 3-RPR (resp. 3-PRR) manipulator^ More recently, studied the RPR- 
2PRR, a simpler planar 3-DOF manipulator that lends itself to algebraic calculus 
||6|. In both papers, the cusp configurations were determined by looking for the 
triple roots of a univariate polynomial. This approach may yield spurious solutions. 
In this paper, the cuspidal configuration are determined directly from the Jacobian of 
the whole set of geometric constraints of the robot, which guaranties that only true 
solutions are obtained. Then, we classify the parameter space of a family of RPR- 
2PRR manipulators according to the number of cuspidal configurations. It is shown 
that these manipulators have either or 16 cuspidal configurations. The proposed 
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method is based on the notion of discriminant varieties and cylindrical algebraic de- 
composition, and resorts to Grobner bases for the solutions of systems of equations. 



2 Mechanism under study 
2.1 Kinematic equations 

A RPR-2PRR parallel manipulator is 
shown in Fig.[T] This manipulator was an- 
alyzed in |6|. It has 1 actuated prismatic 
joint pi, and 2 passive prismatic joints 
P2 and p3. The two revolute joints cen- 
tered in A2 and A3 are actuated while the 
ones centered in Ai, Bi, B2 and B3 are 
passive. The pose of the moving platform 
is described by the position coordinates 
{x,y) of Bi and by the orientation a of 
the moving platform BiBt. The input vari- 
ables (actuated joints values) are defined 
by pi, 02 and 63. The points B\, B2 and 
B3 are aUgned, a= {B\, B2), b= (Bi, B3), 
L2 = (A2,B2)andL3-(A3,B3 




Fig. 1 A RPR-2PRR parallel manipulators 
with (a = 1, = 2, L2 = 2, L3 = 2, X = 1 /2,y = 
1, 9 = 0.2). The actuated joint symbols were 
filled in gray. 



The geometric constraints can be expressed by the following 5 equations 161: 



/2 : X =p2 + L2Cos(62) — ^cos(a) 
/3 :y =L2sin(e2)-i)sin(o;) 



/4 : X —L^ cos(93) — acos(o;) 

fs : y =p3 + L3 sin(93) — asin(a) 



(1) 



Without loss of generality, we fix a = 1 in the rest of the article. 



2.2 An algebraic model 

The singular and cuspidal equations were previously computed using the three fol- 
lowing steps ||6|: 

1 . Reduce the equation system to a polynomial equation depending on the articular 
variables pi , 0i , 62 and one pose variable a: g{pi ,9i,92,(x) =0 (this step is done 
by eliminating variables, either with resultants or with Grobner basis ) 

2. Add the constraint ^ = to define the parallel singularities 

3. Add the constraint ^ = and = to define the cuspidal configurations 

This approach has the advantage of reducing the problem of computing the cusp 
configurations, to the problem of analysing the triple roots of a single polynomial. 
However, this only gives a necessary condition for the manipulator to have cusp 
configurations. In particular it is possible that 3 configurations of the robot coalesce 
in one coordinate but not in the others. 
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Let us come back to the theoretical definitions, using Jacobian matrices to define 
directly the triple roots of the original system of equations in all the input and output 
variables. If P is a list of polynomials and X a list of variables, let Jk{P,X) be the 
union of P = {/^i , . . . , p„,} and of all the kxk minors of the Jacobian matrix of the 
Pi with respect to the X,. 

For the analysis of the RPR-2PRR manipulator, we introduce: 

Y:=[x,y,a,p2,p3\ W:= [fe,L2,i3,Pi,02,03] S -.^ {fu- ■ ■ Js}- 

Using these notations, the parallel singularities of the manipulator are defined by 

{v e M", /:>(¥) =0,^(v) > 0,V/:> £ J5{S,Y)yq G {^,L2,L3,Pi} so that the cuspidal 

configurations are fully characterized by : 

^ = {v e M}\p{y) = 0,^(v) > Qyp £ Js{J5{S,Y),Y)yq G {fo,L2,i3,Pi}} 

3 Main tools from computational algebra 

The algebraic problem to be solved is basically related to the resolution of polyno- 
mial parametric systems. 

More specifically, one needs to solve a system of the following form : 

£ = {v e M",pi(v) = 0, . . . ,p,„(v) = 0,^i(v) > 0, . ..q,{y) > 0} 

where pi, . . . ,p„,,qi, . . . ,qi are polynomials with rational coefficients depending 
on the unknowns X = [Xi,... ,X„] and on the parameters U — [t/i , . . . , C/^/]. 

There are numerous possible ways of solving parametric systems in general. Here 
we focus on the use of Discriminant Varieties (DV, (81) and Cylindrical Algebraic 
Decomposition (CAD, ||3]) for two reasons. It provides a formal decomposition of 
the parameter space through an exactly known algebraic variety (no approximation). 
It has been already successfully used for similar mathematical classes of problems 
(see H). 

To reduce the dimension of the parameter space to three so that it can be dis- 
played, we set L2 = L3. Not that the proposed method can treat the general case 
L2 ^ L3 without any problems. When L2 = L3, the system to solve is ^ with the 
unknowns a,p2,P3, 02, 63] and the parameters [b,L2,pi]. 

3.1 Basic black-boxes 

First experiments are often performed for specific values of the parameters, espe- 
cially singular and/or degenerated cases. Here, we mainly use exact computations, 
namely formal elimination of variables (resultants, Grobner bases) and resolution of 
systems with a finite number of solutions, including univariate polynomials. 

Let us describe the global solver for zero-dimensional systems. It will be used as 
a black box in the general algorithm we describe in the sequel. 

Given any system of equations pi=Q,...,p,„ = of polynomials of Q[A'i , . . . ,X„], 
one first computes a Grobner basis of the ideal < pi,...,p,„ > for any ordering. 

At this stage, one can detect easily if the system has or has not finitely many 
complex solutions. 
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If yes, then compute a so called Rational Univariate Representation or RUR (see 
ifTOl ) of < pi,...,pi„ >, which is, in short, an equivalent system of the form : 
{/(r) = 0,Xi = ^^,...,X„ = ^^}, T being a new variable that is indepen- 
dent of , . . . equipped with a so called separating element (injective on the 
solutions of the system) m G , . . . ,X„] and such that : 

y((pi,. ..,;,,„)) A v{f) ^ v((pi,...,p„,)) 

(xi,...,jc„) I— > ^ = m(xi, . . . ,x„) I—)- ( g'(j3^)' I • • • I g( j3^)' ) 

defines abijection between the (real) roots of the system (denoted by V{pi,... ,Pm)) 
and the (real) roots of the univariate polynomial (denoted by V{f)). 

We then solve the univariate polynomial /, computing so called isolating in- 
tervals for its real roots, say non-overlapping intervals with rational bounds that 
contain a unique real root of / (see ifTTll ). Finally, interval arithmetic is used for 
getting isolating boxes of the real roots of the system (say non overlapping products 
of intervals with rational bounds containing a unique real root of the system), by 
studying the RUR over the isolating intervals of /. In practice, we use the function 
RootFinding[Isolate] from Maple software, which performs exactly the computa- 
tions described above. 

For example, with {b ~ 2, L2 — 2, L3 = 3, pi = 2), the polynomial system ^ 
defining the cuspidal configurations has 16 real solutions. One of these solutions is 
(P2 = 2.30,p3 = 1.86,0 = -2.36,x= 1.79, y = 0.885, 02 = -2.88,03 = -0.999). 
We observe the three coalescing configurations around this root by calculating the 
direct kinematic solutions with 02 = —2.892 and 03 = —1.007. These solutions are 
shown in Figure |2] 




Fig. 2 A cuspidal configuration (left), and the three converging configurations (right). 

3.2 Discriminant varieties 

The above method allows one to study instances of the problem and may be used 
together with a discretization of the parameter space to get a first idea of the com- 
plexity of the general problem to be solved. But the true issue addressed in this paper 
is to find criteria on the parameters that allow classifying the configurations to be 
studied (for example to distinguish the manipulators having cuspidal configurations 
from the others). This leads to a more general problem since one then has to study 
non zero-dimensional, semi-algebraic sets. 

Let pi,. . . ,p,n,qi,- ■ ■ ,qi be polynomials with rational coefficients depending on 
the unknowns Xi,... ,X„ and on the parameters Ui,...,Ud. Let us consider the con- 
structible set : 
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■r-lveC" , /.i(v)=0,...,p„,(v)=0,^i(v)^0,...<?/(v)^0} 

If we assume that ^ is a finite number of points for almost all the parameter 
values, a discriminant variety Vd of is a variety in the parameter space such 
that, over each connected open set '-^ satisfying n Vd = 0, defines an analytic 
covering. In particular, the number of points of '^^ over any point of is constant. 

Let us now consider the following semi-algebraic set : 

^ = {v e Rii,p(v) = 0,^(v) > 0,Vp(v) e J5{J5iS,Y),Y),yq{y) e {/7,L2,i3,Pi}} 

If we assume that has a finite number of solutions over at least one real point 
that does not belong to Vd, then Vd H M'^ can be viewed as a real discriminant variety 
of S^, with the same property : over each connected open set C such that 
^ DVd = V>, '^la defines an analytic covering. In particular, the number of points of 
R over any point of ^ is constant. 

Discriminant varieties can be computed using basic and well known tools from 
computer algebra such as Grobner bases (see fS]) and a full package computing 
such objects in a general framework is available in Maple software through the 
RootFinding[Parametnc] psickage. Figure [5] represents the discriminant variety of 
the cuspidal configurations of the RPR-2PRR manipulator. 

3.3 The complementary of a discriminant variety 

At this stage, we know, by construction, that over any simply connected open set 
that does not intersect the discriminant variety (so-called regions), the system has a 
constant number of (real) roots. 

The goal of this part is now to provide a description of the regions for which the 
number of solutions of the system at hand is constant. For that, we compute an open 
CAD (El E]). 

Let CQ[Ui, . . . , Uii] be a set of polynomials. For ; = c/ — 1 . . . 0, we introduce 
a set of polynomials C Q[t/i , . . . , Uj^i] defined by a backward recursion; 

• : the polynomials defining the discriminant variety 

• : {Discriminant(/9,f/,),LeadingCoefficient(;?,t/,), Resultant(p,^,f/,), 

We can associate to each an algebraic variety of dimension at most ; — 1 : V,- = 
V{Y[peiPiP)- Figure[3]and|4]represent respectively V3 and V2 for the manipulator at 
hand. 

The Vi are used to define recursively a finite union of simply connected open 
subsets of R' of dimension /: U^l^j^- jt such that V,- n '^■jt = 0, and one point Uij^ 
with rational coordinates in each ,t. 

In order to define the '^n^, we introduce the following notations. If p is a univari- 
ate polynomial with n real roots: 

{-ooif / <0 
the real roots of p if 1 < / < n 
+00 if / > n 
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Fig. 3 Discriminant Variety of the cuspidal Fig. 4 V2 for the cuspidal configurations of 

configurations the manipulator. 

Moreover, if /? is a «-variate polynomial, and v is a « — 1-uplet, then denotes 
the univariate polynomial where the first « — 1 variables have been replaced by v. 
Roughly speaking, the recursive process defining the "Wi k is the following: 

• For ; = 1, let /?! = Y[pe.iPi P- Taking '^i j^ =]Root {p, k); Root {p,k+l)[ for k from 
to n where n is the number of real roots of pi, one gets a partition of R that fits 
the above definition. Moreover, one can chose arbitrarily one rational point mj 
in each "^i jt. 

• Then, let pi = IlpGiifi P- The regions ^jj^ and the points m, are of the form: 

%k = {(vi,...,Vi_i,v,) |v:= (vi,...,v,_i) e ^-ij, 

ViG]RootipJ,l),RootipJJ + l)[} 
I a a a\ -^u f (^b A-l) = "i-l./ 

=(A,...,A-i,A), (^,.e]/^oor(p;."-'•^/),Lor(p;."-'■^/+l)[ 

where y',/ are fixed integer. 

For our example, we get for 773 a trivariate polynomial of degree 33, for p2 a 
bivariate polynomial of degree 113, and for p\ a univariate polynomial of degree 
59. The zero-dimensional solver then provides the positive real roots of p\ (Table 
[Til, from which we easily deduce the open intervals uw^. We then use the zero- 
dimensional solver to solve every pi(u\±-,Vi) and deduce all the tests points of 
the in each cells of Figure]?] Finally we use the zero-dimensional solver to 
solve every P3{u2,kiU3) and deduce all the tests points of the '^3 J./, describing so 
the complementary of the discriminant variety. 



Table 1 Numerical values of the positive roots of pi 



b 


b, 


bi 


bi 


b4 


bs 


be 


bi 


b» 


bg 




bu 


bi2 




0.0 


0.533 


0.564 


0.617 


0.656 


0.707 


1.0 


1.41 


1.52 


1.62 


1.77 


1.88 



3.4 Discussing the number of solutions of the parametric system. 

At this stage, we have a full description of the complementary of the discriminant 
variety of the system to be solved : a recursive process for the construction of each 
cell (t and a test point (with rational coordinates) in each of these cells. By def- 
inition of the discriminant variety, we know that the system has a constant finite 
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number of solutions over each of these cells and computing this number for each 
cell is the only remaining step. This can be done simply by solving all the systems 
■^\v=uj r ^ ~ 1 • • • "f/ using the zero-dimensional solver. 

For our example, the process described in [33] returns 344 cells of dimension 3 
( '^31, ...,'^3,344 ). We solve the system S'' for each of the 344 associated sample 
points, and we get always either or 16 solutions. By selecting only the cells where 
the manipulator has 16 cuspidal configurations, we obtain the 58 cells shown in Fig- 
ure |5] Table |2] provides the different formula bounding the three dimensional cells 
^^3,1, "^3344 and Table [3]represents the 58 cells of Figure|5] where the manipula- 
tor has 16 cuspidal configurations. 



Table 2 Formula describing the boundaries of the cells in Table[3] 



fcl =0.Z>2 =Rool(86^ - llfc'^ + 6/?2 -1,2) 

^3 = Root{'ib^+b^ -3b^ -l,2},b4=Ro<j\{b^ +3b^ +ib^ -i.2) 

b^=Root{-2b^+b^+3b^-l,l).bf^ = 1/V2,i>7 = hbg = ^2) 

69 = Root(262 +b^-3b^-\.2),biQ= Root(Z7^ - ft^ - 3 // -3 fo- - 1 , 2) 

b^i =Root(-4fc^+i^ + 3fc2-L2) 

il2 =Rool(fo^-6fo'*+ 11 6^-8,2), i>j3 =00 


i2| 1*)= RiJ0l(^(6- + iy'L^-3(6-+3/j+l) (ii- - 36+ l) (fc- 1)2 (6+ 1)2^ + 

3(^b^ + lj {b-\f (6+l)''t2_(6_i)6(j^[)6 2^ 
i22(i>) = 1 -6^^23(6) = = (1 -6-)/6.t25 (*) = •« 
Z.2^ (6) = b^/Vl-b^, L27 (6) - 1 / Vb^-l. L2g (6) = i^/v/fcJ - 1 
hg(l>)= {l>^-n/l',L2^g{b) = b'-l 


Plj (i.L2) = R<KJl(-pfi>^ +36* (l^S^ + I - t?)pf -36^ (^-7 L^b- + 7 1^+I^+I^b^ -lijb- + lj p J + + 1 -L^)"' ,2) 
Pl2(^^4) = ^'>°t(pf + (-36'^+3L^/)--3Z.^^ pj+ (llL^b^ +3 L^b* -f,L^b~ +2 b^ -21 L^b"^ +3 L^'jp^ + (^L^b^ -b-^ -L^y .2) 
Pl^(i.Z,2) = i-,pi^(6,L2) = l/6 



Table 3 Cells of M where the manipulator has cuspidal configurations. 



I61621 


ll''2, i22 l-IPl 1 PI2 D' ll'-22 '■23 UPI3 Pl, D- (li23 UPlj Ph !)• (1% % UPI3 PI4 1) 


1*2 '3 [ 


(lt2,i2jl-lPl,Pl2l)'(l'-26'-22HPl,Pl2l)-(li22'-2,Upl3Pl2l)-(|i23i24Upl,Pl2l)'(l'-24i-25Upl,Pl4l) 


Ihhl 


(|t2ii26l-IPliPl2U^(l'-26'-22MPl,Pl2U-(|i22''24UPl3Pl2lM|i24t2,l-IPl,Pl4l)-(li2,'-25MPl,Pl4l) 




(lhj'-2^l-lPllPl20:U\h2l-lPliPl2i'>-fi'<i2'^-4^-'l'h''hJ^''f''l4''l}^^^ 


ihi'el 


QLl, i2,l-lPliPl2l)-(li22'-26UPl3Pl2l).(|i26''24UPl3Pl2l)-(|i24i2i UPl, PI4 I)- 0^2, ^25 UPl, PI4 1) 


Hhl 


a'^2, i22 l-IPl 1 Pi, 11- ll''22 % I'lPlj PI2 % UPI3 PI4 !)■ tihf, % UPI3 PI4 !)■ (l'-23 ^25 l-IPlj PI4 D 


Ihhl 


'ih, i29l-IPl2Plill'll''29'-2ioHPl4Plil''(l'-2io'-27HPl4Pl3l)-(li27'-28l-IPl4Pl3l)-'!i28i25l-lPl4Pl3l) 


]68 69 [ 


QLj, i29l-lPl,Pl,[)-(li-29 i27[.lPl4Pl, [)■ (li27 ^2,0 [■IPI4 Pi 1 D' (fe^ i2g l-IPl4 Pi, !>■ (.^j 4, .■ IPIj PI, [I 


I69610I 


(|L2jL2,|.|Pi,Pi,l),l|i27'-29l-IPl2PliU-(|i29i2,oHPl4PliU'll'-2,o'-28l-lPl4Pl,U-ll'-2gi2,l-IPl4Pl3l) 


I'loi'lll 


(|t2, i27 l-IPl, Pi 1 1)^ (1*^27 '■29 UPI2 Pi 1 1). (|i29 i2,o UPI4 Pi 1 »■ (l'-2io l-IPl* Pi, O-Uhg % l-IPlj Pi, D 


]^11 ''12[ 


(]L2, i27 l-lPl, Pi 1 1). (1^27 ^29 UPb Pi 1 D- (|i2, i2g UPI4 Pi , I). (|i28 ^210 HPI4 Pi , I- Oi2io '■25 l-IPlj PI3 D 


]^12''13[ 


(li-2i i29 l-lPl, Pi 1 [)' Ili29 i-2g [■IP^ Pi , [)■ (li28 ^2,,, l-IPl4 Pi , [)■ (fejo ^25 l-IPl4 PI3 [) 



4 Conclusion 

We have proposed a general method to describe rigorously the design parameters 
for which a manipulator has cuspidal configurations. This method can be applied 
directly to other mechanisms, such as the ones studied in 141 [121 for example. The 
tools used to perform the computations were implemented in a Maple library called 
Siropc0. For 3D illustration purposes, we have detailed the main computations to 
be performed with manipulators satisfying L2 = ^3. However, the proposed method 
allows one directly to solve the general case (L2 7^ ^^3) by computing a discrimi- 
nant variety of the system with 4 parameters b,L2,Lj,,pi, and by decomposing 
with a CAD adapted to the discriminant variety. This description generalizes and 



http://www.irccyn.ec-nantes.fr/ moroz/siropa/doc 



8 



G. Moroz'-2, D. Chablat', P. Wenger', F. Rouiller^ 



completes the analyse done in I6l. There is still some limitations though. In particu- 
lar, when the system that defines the cuspidal configurations has no solution, it may 
mean that there exists a manipulator with no cuspidal configurations, but it may also 
mean that no manipulator can be assembled with these design parameters. Thus it 
is essential to be able to describe precisely the set of design parameter values for 
which a manipulator can be assembled. 





Fig. 5 The cells of M? where the manipulator admits cuspidal configurations, front view (left) and 
back view (right). 
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